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Abstract 

We interpret the Lorentz force equation as a geodesic equation associated with a non- 
linear connection. Using a geometric averaging procedure, we prove that for narrow and 
smooth one-particle distribution functions whose supports are invariant under the flow of the 
Lorentz equation, a bunch of charged point particles can be described by a charged cold fluid 
model in the ultra-relativistic regime. The method used to prove this result does not require 
additional hypotheses on the higher moments of the distribution. This is accomplished by 
estimating the expressions that include the differential operators appearing in the charged 
cold fluid model equation. Under the specified conditions of narrowness and ultra-relativistic 
dynamics, it turns out that these differential expressions are close to zero, justifying the use 
of the charged cold fluid model. The method presented in the work can also be applied to 
justify the use of warm plasmas and other models. 



1 Introduction 



Despite limitations concerning the mathematical description of the discrete nature of the parti- 
cles comprising a plasma, modeling the dynamics of relativistic non-neutral plasmas and charged 
particles beams by fluid models is a common place. The relative simplicity of these models (com- 
pared with the corresponding kinetic models) is a partial justification of their use. We propose 
in this work another justification. We will concentrate in the charged cold fluid model. However, 
we should notice that the same philosophy is applicable to more sophisticated models. 

In high intensity beam accelerator machines, each bunch in a beam contains a large number 
of identical particles contained in a small phase-space region. In such conditions, a number 
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of the order of 10 — 10 charged particles move together under the action of both external 
and internal electromagnetic fields. Often in modern applications, such beams of particles are 
ultr a-relativistic . 

One is interested in modeling these physical systems in such a way that: 

1. The model for a bunch of particles must be simple, in order to be useful in numerical 
simulations of beam dynamics and for analytical treatment, 

2. It allows for stability analysis and a qualitative understanding of the dynamical behavior 
of the system. Three dimensional numerical simulations can be also desirable. 

The standard approach has been to use fluid models as an approximation to a kinetic model. 
These derivations of fluid models from kinetic models can be found in [1-4] and references there. 
They are based on some assumptions, usually in the form of equations of state for the fluid or 
assumptions on the higher moments of the distribution function. This constraint is necessary 
to close the hierarchy of moments of the distribution function and to have a sufficient number 
of differential relations to determine the remaining moments. This is a general feature of all 
the derivations of fluid models from Kinetic Theory (see for instance [1-4] and references there) , 
where a truncation scheme is required for the fluid model to be predictive. 

We present in this work a new justification of the charged cold fluid model from the framework 
of kinetic theory. The novelty of the new approach is that it uses natural hypotheses suitable 
for particle accelerator machines and exploit only the mathematical structure of the classical 
electrodynamics of point particles interacting with external electromagnetic fields. In particu- 
lar, we estimate the covariant derivative of the mean velocity calculated with the one-particle 
distribution function. This is given as an asymptotic formula in terms of the time of evolution, 
diameter of the distribution and energy of the beam. The fluid model is described by only one 
dynamical variable, the normalized mean velocity field. The variance and the heat fiow tensor 
are not necessarily zero but are finite. In our treatment both the fluid energy tensor and the 
fluxes tensor are non-dynamical. This is why we have considered here the charged cold fluid 
model. However, the approach is consistent, since our aim is not to give an equation for the 
mean velocity field, but evaluate how much certain differential expressions (formally equivalent 
to the charged cold fluid model equations) differ from zero. One can stipulate the validity of the 
model from the estimate of the corresponding differential expressions. On the other hand, in the 
models presented for instance in [1-4], the variance and the covariant heat flow arc dynamical 
and a system of partial differential equations is used to determine the dynamics of these fields. 
However, assumptions must be done to close the system of fields equations and the validity of 
those assumptions is unclear. 

The method used to obtain these bounds is the following: 

1. Firstly, we re-write the Lorentz force equation as an auto-parallel condition of a non-affine 
linear connection. 

2. We use an averaging procedure described in [5] to average this connection. 

3. The resulting averaged connection is an affine connection on the manifold M. 

4. It happens that under the same assumptions as used for the particle dynamics [6], the 
corresponding solutions of the Vlasov equation / and the averaged Vlasov equation / are 
similar. This comparison result is based on the comparison results of the point dynamics. 
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5. Given the distributions / and /, one has the corresponding mean velocity vector fields. 
One can prove under the same hypotheses that these mean velocity vector fields are similar, 
which means that the difference between them is controlled by powers of a small parameter. 

6. Finally, we show here that the auto-parallel condition of velocity field of the averaged 
Vlasov equation associated with the averaged dynamics is controlled by the diameter of 
the distribution /. Together with the above point, this result provides estimates for the 
auto-parallel condition of the mean velocity field of the solution of the Vlasov equation. 

Therefore, the methods presented here and the usual derivations of the fluid models con- 
tained in [1-4] are different. The standard approaches assume an asymptotic expansions of the 
differential equations for the moments, in terms of a perturbation parameter which is similar to 
the diameter a of the distribution function. These asymptotic expansions are realized at low 
orders as a truncation scheme in hierarchy of moments. Then, they discuss systems of partial 
differential equations such that they arc self-contained and consistent with physical constraints 
and with the asymptotic expansions. On the other hand, our approach is based on the structure 
of the Lorentz force equation of a charged point particle, which lies at the basis of the kinetic 
models. After being written in a geometric way, the Lorentz force equation is substituted by the 
averaged Lorentz force equation. The key points are that the averaged Lorentz force equation 
depends only on the first, second and third moments of the distribution function and that is 
a simpler equation than the original Lorentz force equation. Then, under some smoothness 
assumptions on the distribution function, we can place bounds on the differential expression of 
interest. 

There are some assumptions that we have used in the present work. For instance, we assume 
that the distribution functions have compact support. This excludes important examples like 
Gaussian distributions, which are the distributions in equilibrium. However, one can consider 
truncated Gaussian distribution as approximations, when the truncation is performed smoothly, 
using a bumpy functions techniques. 

Also we have assumed that the distributions functions are smooth (at least of class C^). 
Although we do not currently have a proof that we can extent our results to bigger functional 
spaces for the distribution functions, since the main results are written in terms of Sobolev 
norms, it is conjectured that they can be extended to Sobolev spaces. 

Notation. 

Let (x, U) be a local coordinate system on M , where U C M is an open sub-set and 
X : U ^ R" a local coordinate system. An arbitrary tangent vector at the point p G U 
is of the form Xp := X = X^-^\p. The local coordinates associated to the tangent vector 
Xp € Tj^M C TM are {x^,y^). M will be n-dimensional and equipped with a metric of 
signature —). We call M space-time manifold. In some calculations it will be useful 

to consider the particular case when the metric is flat. Given a 1- form (w" means duality 
defined by the semi-Riemannian metric 77. 

There are certain sub-bundles of the tangent bundle that are relevant in our treatment: 

1. The slit tangent bundle N := Uxgtm iv ^ T^jM, rnj^x) y'^y^ > }, 

2. The unit hyperboloid bundle S := U^eTM iv ^ T^M, r]ij{x) y'y^ = 1 } 

3. The null bundle NC := Ugtm {2/ ^ T^M, Vij{x)fy^ = 0}. 

The particularization to each point x G M are and E-^ respectively. 
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We are dealing with one-particle distribution function f(x,y). We assume that fx '■= fix,) 
has compact support on the unit hyperboloid bundle X). The diameter of the distribution fx is 
ax ■= sup{dfj{yi,y2) 1 2/1,2/2 S supp{fx)}- Then, we define a := sup{ax, x G M}. We define the 
energy function of a distribution / to be the real function: 

E:M — ^ R 

X E{x) := inf{y°, y G supp{fx)}, 

where is the 0-component of a tangent vector of a possible trajectory of a point particle, 

measured in the laboratory coordinate frame. 

We introduce from reference [6] our notion of semi-Randers space. 

Definition 1.1 A semi-Randers space consists of a triplet (M,?7, [A]), where M. is a space-time 
manifold, rj is a semi-Riemannian metric continuous on M and it is smooth on TM \ NC and 
the class of locally smooth 1- forms A G [A] defined such that F = dA for any A €F. 



2 Lorentz Force Equation and Averaged Lorentz Force Equation 

We proposed in [6] a geometric description of the dynamics of a charged point particle interacting 
with an external electromagnetic field. The relevant data was extracted from a Lorentzian metric 
r] with signature (+,—,—,—) and from the Lorentz force equation, that in an arbitrary local 
coordinate system reads [7,8]: 



d^cT* „ . da^ da^ ,-,/,,x da^ I .da da, 

where a : I — > M is a solution, a curve parameterized by the proper-time r associated with 
r/, jk are the coefficients of the Levi-Civita connection ''V of r] and dA = F is the exterior 
derivative of the 1-form A. 

We viewed the system of differential equations (1.1) as the auto-parallel condition of a con- 
nection that we called the Lorentz connection ^D. In reference [6] were defined several related 
connections: the non-linear connection [9,10], the linear connection in the sense of Koszul [11] 
on a suitable sub-bundle of TTM and an example of a connection in the pull-back bundle 
7r*TM [7]. These connections are related. For instance, the non- linear connection determines 
the linear connection in the sense of Koszul. The non-linear connection (with some additional 
assumptions on the torsion-type tensor) also determines the linear connection on the pull-back 
bundle. In this work we consider mainly Koszul connections (which we will denote by and 
< >) as determined in reference [6], since this allows us to introduce covariant derivatives 
and auto-parallel equations. 

Let us denote hy r]{Z,Y) := r]ij{x)Z^ , 

Definition 2.1 For each tangent vector y € T-jM with r]{y,y) > 0, there are defined the fol- 
lowing functions: 

^r%fe(x,2/) = ^r,fc + -^^=(F%(x)y'"r?^fe + Ffe(x)2/"^r?^,)+ 

-\-F^rn{x) — , (riik 7 TVisVkiy^y''), (2-2) 
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where r]{y,y) is an abbreviation of r]ij{x)y^y^ , ^T^ jk, ii,j,k = 0, 1,2, ...,n) are the connection 
coefficients of the Levi-Civita connection ^'V in a local frame, Fij := diAj — djAi and F* j = 

Proposition 2.2 The Lorentz force equation can be written as the auto-parallel condition 

""D-i = 0, 

where x : I — > M is a time-like curve parameterized with respect to the proper time of the 
Lorentzian metric r], x is the horizontal lift on N and is the Koszul linear connection deter- 
mined with connection coefficients given by equation (2.2). 

Proof: One can check this fact by direct computation (for instance see [6]). □ 

We win define averaging procedure for several geometric objects. These averaging operations 
correspond to fiber integrations of sub-bundles of TM. (in particular we are interested in objects 
living on the unit hyperboloid sub-bundle). Usually the measure is given as f{x,y)dvol{x,y), 
where f{x, y) is the one-particle probability distribution function [12]. The volume form dvol{x, y) 
is induced by the Lorentzian metric rj: 

dvol{x, y) = \/ —detr] dy^ A • • ■dy^~^. 

One can prove the following 

Proposition 2.3 The averaged connection of the Lorentz connection defined on the pull- 
back bundle 7r*TM ^ N is an affine, symmetric connection on M. The connection coefficients 
are given by the formula: 

1 1 

Each of the integrations is equal to the y-integration along the fiber: 

9/(S^) := / f{x,y)dvol{x,y), < y' >:= J" / y' f{x,y) dvol{x,y), 

< y'y^ >:= j [ yV' f{x, y) dvol{x, y), 

VOl[2^x) JEa, 

< yV^'/ >:= ,L X / yVy'' fix, y) dvol{x, y). 
voliT,^) Js^ 

Definition 2.4 There is a Koszul connection defined on M such that given sections X,Y E 
TTM is defined by: 

< ^L>>xi^:= (X^^+< ■^r%fc>(x)X*y^)9„ i,j,k = 0,--;n-l. (2.4) 
We call this connection the averaged Lorentz connection. 

Proposition 2.5 Let M be the space-time manifold and A a 1-form on M. Assume the existence 
of a non-negative function f : S — ^ R with compact support. Then 



vol 
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1. The connection < > is an affine, symmetric connection on M. Therefore, for any 
point X G M, there is a normal coordinate system such that the averaged connection 
coefficients are zero. 

2. < > is determined by the first, second and third moments of the distribution function 

f{x,y)- 

Later wc will obtain boTinds on the value of certain covariant derivatives. These bounds are 
obtained using a norm, defined in a particular local frame. This local frame is associated with 
the bulk motion of the fluid. In particular the local frame is determined using the mean velocity 
field <f>i,. 

In a Lorentzian manifold (M, rj) where there is defined on M a time-like vector field U nor- 
malized such that r]{U, U) = 1, one can define the Riemannian metric fj : 



fj{X, Y) := -n{X, Y) + 2n{X, UUY, U). (2.5) 

We choose the following locally smooth vector field U in the definition of the Riemannian metric 
(2.5): for r}ij{x) < > {x) < > {x) > 0, U{x) is given by 

U(x) = <y>i^) (2 6) 

^/Vijix) <t> {x)<y^ > (xY 

This vector field is not continuous, since there is a discontinuity in the boundary d(7r(supp(f)). 
Also let us note that for at each x and for functions / with convex supports on the connected 
hyperboloid, the mean velocity field is time-like. From now, we will assume that the support of 
the distribution function is connected. 

Using f] there is a Riemannian metric on each vector space Ta;M defined as fiij{x)dy^ dyK 
It induces a distance function dfj on the manifold T-^M and related operator norm [13]; given a 
continuous operator Ox : Tj;M — > T^jM, its operator norm is defined by 

IIO.II, := sup{ ^^f^(x), y G T,M \ {0} |. 

^ \\y\\ri ^ 

Let us restrict to the case where the Lorentzian metric is the Minkowski metric in dimension 
n. Let us denote by 7(t) the gamma factor of the Lorentz transformation from the local frame 
defined by the vector field U to the laboratory frame, at some instance defined by the local time 
t, the coordinate time defined by the laboratory frame. Denote by O'^(t) = 'y^{t)— < y >^ (t) 
and 9'^{t) =< y (t) — y^{t). Here y(t) is the velocity tangent vector field along a solution 
of the Lorentz force equation and y{t) is spatial component of the tangent vector field along a 
solution of the averaged equation, with both solutions having the same initial conditions. The 
maximal values of this quantities on the compact space-time manifold are denoted by 6^ and 9^. 

The basic result from [6] that we use here is 

Theorem 2.6 Let (M.,ri,[A]) be a semi-Randers space and rj the Minkowski metric. Let us 
assume that: 

1. The auto-parallel curves of unit velocity of the connections and < > are defined 
for the time t, the time coordinate measured in the laboratory frame. 

2. The dynamics occurs in the ultra-relativistic limit, E{x) » 1 for all x G M. 

3. The distribution function is narrow in the sense that a « 1 for all x eM. 
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4- It holds the following inequality holds: 

-^1 < 1, 

5. The support of the distribution function f is invariant under the flow of the Lorentz force 
equation. 

6. The following adiabatic condition holds: ^^^^^ << 1. 

Then, for the same arbitrary initial condition {x{Q),x(0)), the solutions of the equations 

^VxX = 0, < >j X = 

are such that: 

\\x{t) - x{t)\\ < 2{C{x)\\F\\{x) +Ci{x){l + B2ix)a))a^E-\x)t^ (2.7) 
where the functions C{x), Ci{x) and Bi{x) are bounded by constants of order 1. 

In a similar way, one can compare the velocity tangent fields along the corresponding geodesies 
[6]: 

Theorem 2.7 Under the same hypothesis as in theorem 2.6, the difference between the tangent 
vectors is given by 

\\3:{t)-x{t)\\ < {K{x)\\F\\{x) + KI{1 + D2{x)a))a'^ E-^ t. (2.8) 

with Ki and D2{x) functions bounded by constants whose values are unity within an order of 
magnitude. 



3 Comparison of the Solutions of the Vlasov and Averaged 
Vlasov Equations 

In this section we estimate the difference between the solutions of the Liouville equations asso- 
ciated with the averaged Lorentz connection and the original Lorentz connection. 

Given a non-linear connection characterized by the second order vector field x ^ TTM, the 
associated Liouville equation is %(/) = 0. 

Examples. 

1. From the connection Lorentz connection ■^V*^(x,y) [6], one can recover the spray coeffi- 
cients ^G{x, y), using the homogeneous properties on y and Euler's theorem. In particular, 
the spray coefficients are: 

^G*(a;,y)= ^r%fc(a;,2/)y^/= hv ^„ + -j2=={F^ ,{x)y"'rjmk + F\{x)y"^Vmj)+ 

= hr^j, + (F%(x)y-r?^fc + F*fc(x)y™7?^,))y^/ = 

V 2y/rj{y,y) ^ 
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Prom here, one defines the vector field Lx' 

2. A similar procedure applies to the averaged Lorentz Vlasov vector field. In this case, 
however, we do not have the simplifications above, except for rj(y,y) = 1. Therefore, the 
spray coefficients are: 

< > {x,y) = '^rjk{x,y)y^y'' = jk+ < -^^={F^ j{x)y"'Vmk+F\{x)y"'vmj)+ 

^ Wv{y,y) 
+F' m{x)—^=={vjk - -f — TVjsVkw'y^)) > y^y^ = 

+F'mix)—i7]jk - rjjsmiy'y^)) > y^y^ = 

= (T^-fe + \{F'j{x) <y^> rimk + F\{x) < y'" > Vmj)+ 

+F'^{x){r]jk <^> -Vjsrjki < ^y'y' >)) y^y''- 



The averaged Vlasov vector field can be written in a similar way as before, 

d . ^ 1 



< X >= y 



dx'^ 



("r^-fey^/ + -(F%(x) <y^> vmk + F\{x) < y^ > Vmj)+ 



+F'm{x){vjk < ^ > -VjsVki < ^y"y^ >)) y^y''- 

Proposition 3.1 Let f and f be solutions of the Vlasov equation ^xif) = ^^e averaged 
Vlasov equation <^ x > (/) = 0? where ^x ^.f^d < ^X > C'f^ ^he spray vector fields obtained 
from the non-linear connections and < >. Let us assume: 

1. The same hypotheses as those in theorem 2.6, 

2. supp{f) C supp{f), 

3. supp{f) is a sub-manifold o/TM. 

Then for the solutions of the Vlasov and averaged Vlasov 's equation one has the relation: 

\f{t, x{t),x{t)) - fit, x{t),xm < {C{x)\\F\\f^Ci{x){l + B2{x)a)) E'^ t^ + 

+ {k{x)\\F\\r,{x) K^{1 + D2{x)a))a'^ E-'^ t. (3.1) 
for some functions C{x{t)) K[x{t)) along the geodesic of the Lorentz connection. 
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Proof: / and / are solutions of the corresponding Vlasov and averaged Vlasov equations re- 
spectively. Therefore, / and / are constant along the corresponding auto-parallel curves: 

\f = j^f{x{tlx{t)) =0, <\>f = jjix{t),m = 0. 

For the same initial conditions, the geodesic curves corresponding to the connections ■'^V and 
arc nearby curves at the instant t in the way described by theorem 2.6. 

Let us introduce the family of interpolating connections, 

^Ve :=(l-e)^V + e< W>. 

Each of them has an associated spray vector field x^, e G [0,1]. Therefore, let us consider a 
solution fe{x, y) be the solution of the following Liouville equation: 

^XJ. = 0. 

Since the dependence on (e,x, y) of the vector field ^Xe is C^, the solutions of the differential 
equation are Lipschitz on the parameter e. We can see this fact in the following way. The 
Liouville equation can be written as: 

^Xje = ^^/(X,(s),i/e(s))=0, 

where {x^{s),ye{s)) is an integral curve of the vector field ^Xe restricted to the unit hyperboloid 
bundle and such that it is parameterized by the coordinate time t. 

One can use standard results from the theory of ordinary differential equations to study 
smoothness of the solutions of the above equation. In particular, the connection coefficients for 
the interpolating connection are 

( ^r,)^ := (1 - e) jk+e< "^V jk > ■ 

Prom the formula (2.3) for the coefficients ^P* jjt one can check that (^Fg)* are smooth func- 
tions in the open set of time-like vectors y and on the parameter e. 

We will give an upper bound on the difference \f{t,x{t),x{t)) — f{t,x{t),x{t))\. In order to 
achieve this, standard results on the smoothness of the solution of differential equations are 
used (sec the chapter 1 of [14]). In particular we use that for each (e, x{s),y{.s)), there is an 
open neighborhood Ue of [0, 1] x supp{f ) containing (e, x{t),y{t)) such that the solutions of the 
differential equations are Lipschitz in Ug: 

\r{t,xe{t),xe{t)) - f{t,xe{t),x-,{t))\ < ci{e, x{t) , ^t))6{{e, x{t) , 3i{t))) + 

+C2{e,x{t),x{t)) \\Xe{t) - Xe{t)\\fj + C3{e,x{t),x{t)) \\Xe{t) - Xe{t)\\fj. 

Ci(e, x{t),x{t)) are constants which depend on the open neighborhood U^; 5{{e, x{t), x{t))) is the 
diameter on the e component where we are applying the Lipschitz condition. 

One can always choose a refinement of an open cover of [0, 1] x supp{f) such that both the 
Lipschitz condition, theorem (2.6) and theorem (2.7) can be applied simultaneously. Since [0, 1] 
and supp{f) are compact, we can consider a finite open covering of [0, 1] x supp{f). Then, using 
the above local bound one obtains a global bound of the form: 

\fit, Xit),xit)) - f{t,x{t),iit))\ < CI + C2 Mt) - Xit)\\f^ + C3 \\xit) - i{t)\\f^. 
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The constants q are finite (by definition of Liptschitz and by compactness of the interval [0, 1]). 
The functions / and / are constant along the respective geodesies when parameterized by the 
corresponding proper time. Therefore, 

\f{t,x{t),x{t)) - }{t,x{t),im = |/(0,a;(0),x(0)) - /(0, 5(0), 5(0))|. 

Let us assume the same initial conditions x(0) = x(0) and x(0) = a:(0). Since the difference 
\J{t,x{t),x{t)) — f{t,x{t),x{t))\ is a smooth function of \\x{t) — x{t)\\f^ and \\x{t) — x{t)\\f^, one 
has that 

< ci < n • Ki\\x{t) - x{t)\\f^ + Ki\\x{t) - x{t)\\f^ 
for some constants Ki. Then we have that: 

\f{t,x{t),x{t)) - f{t,x{t),x{t))\ < 

< \fit,xit),x{t)) - f{t,x{t),im + \f{t,x{t),x{t)) - fit,xit),i{t))\. 

The first term is bounded by ci, which is bounded by n- Ki\\x{t) — x{t)\\fj + -Ki \\x{t)—x{t) \\f^. The 
second term can be developed in Taylor series on the differences ||a;(t) — x(i)||jj and ||i;(t) — x(t)||jj, 
since / is smooth. Therefore: 

\f{t,x{t),x{t))-~f{t,x{t),xm < {C{x)\\F\\Cl{l + B2{x)a))a^E-H^ + 
+ {k{x)\\F\\{x) Kl{l + D2{x)a))a^ E'H. 

□ 



4 The Charged Cold Fluid Model from the Averaged Vlasov 
Model 

Definition 4.1 Given a semi-Randers space (M, r/, [A]), the averaged Maxwell- Vlasov model is 
defined by the dynamical variables (F, /) determined by the coupled system of equations: 

<^X>f = 0, ^= - dvol{x,y)yf{x,y), p= ^ dvol{x,y)f{x,y). (4.1) 

In this case, the remaining dynamical variable f{x, y) is determined by an integral-differential 
equation: 

<''X>/ = 0, (4.2) 

where is the Liouville vector field of the averaged Lorentz dynamics associated with the 
external electromagnetic field F. 

Since we will use the results of section 3 and reference [6], the distribution function / is at 
least of type C^. Since the support is compact, there are several Sobolev norms [13] which are 
defined. We will write our results in terms of these norms. 

There is one question which deserves to be clarified. Firstly, we have assumed that the 
supp{fx) is a smooth manifold. In principle, this is not necessarily the case. However, we have 
restricted to the case where supp{fx) C supp{fx)- This can be achieved in the following way. 

Proposition 4.2 Let < ^x> fi^jU) = ^xfi^^u) = be such that 

1. The domain of definition of the vector < > is an open sub-manifold ofS. 
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2. supp{fx) C supp{fx) for all x eM, 

3. supp{fx) and d{supp{fx)) are sub-manifolds of S-^ for all x G M. 

Then one can reduce supp{fx) — > supp{fx) for all x G M. 

Proof: Let us consider the product of functions fg, where < > / = and the function 
5f is a bumpy function defined on the domain of the vector field < >■ Since both, the 
domain of definition of the ri(< ^x >) supp{fx) are sub-manifolds, this function exists [15]. 
Therefore, 

^(a:;, y) = 0, (x, y) G supp{f ) \ supp{fx), g{x, y) = 0, {x, y) G dsupp{fx) 
and all the derivatives are zero on dsupp{fx)- Then, one can perform the following calculation: 

< > Cfg) = g<\>f +f<''V>g = 0. 
Therefore, one can always restrict the solutions of < > / = in such a way that supp{fx) = 

SUppifx). □ 

Also note that using C°° bumpy functions we can obtain a distribution function / which is 
bounded by 1: since supp{f) is compact, we can consider the positive function -j-Jj — /, with 

< > / = 0- Then, < ^x > TTiTi — f — ^ the new solution is bounded by 1 (here 

ll/lloo 

ll/lloo := max{f{x,y), y G supp{f)}). 

In the calculations below we assume that supp{ fx) is orientable, with a volume form dvol{x, y). 
If this is not the case, it follows from propositAon 4-2 that we can approximate {fx, supp(fx)) by 
{fx, supp{fx)) and then tisc dvol{x,y) induced from the volume form on the unit hyperboloid 

Hx when wc calculate the averages. 

Using equation (3.1), it follows that the error induced by the substitution / — > f is of order 
and due to proposition 4-2 not harmful reducing the support such that supp{fx) — > supp{fx)- 
Therefore, in the following calculations, when it is convenient we can use dvol{x, y) as a measure 

and substitute supp{fx) by supp{fx)- However we will write the calculation assuming that 
supp{f) is a manifold with a volume form dvol{x,y), since this is a more general statement. 

Let F be a closed differential 2-form defining the Liouville vector field ^x- Let us consider the 
Sobolev spaces {y\^^'^{supp{fx)), \\ ■ ||i,i) and {yV'^'^{supp{fx)), \\ ■ ||o,2) [13]. Recall that the space 
of smooth functions F{supp{fx)) is a sub-space of both of those Sobolev spaces (for compact 
supp{fx))- 

Let us denote by 5{x, y) =< y > {x) — y. 

Theorem 4.3 Let M be an n- dimensional space-time manifold and < > the vector field 
associated with the averaged Lorentz force equation. Assume that: 

1. The distribution function is such that fx, djf G J^{supp{fx)) C yVi,i{supp{fx)) , 

2. The function 6{x, y),djd{x, y) G WQ,2{supp{fx))- 

Then 

1 

II < ^D>yV{x)\\,< ^^f'^PPf^^ \\dolo9{6'x)\\o,2) -Wflhi-a^ +0{a\ (4.3) 

VOl{supp[fx)) ^ 
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where 5x{-) ■= 6{x, •) and 

V\x) =< f >f (x) := ^ / _ dvol{x,y) f{x,y)y\ 

^ LuppCf). y)dvol{x, y) Jsuppif), 

The volumes are the following 

vol{supp{fx)) ■■= vol{T,x); volE{supp{fx)) ■= _ dvol{x,y)\ 

J supp(fx) 

the derivative in equation (4.3) refers to the local frame such that the vector U = {Uo,o) 

Proof: Because the averaged Lorentz connection is an affine connection on M, given a point 

X € M, there is a coordinate system where the connection coefficients are zero at that point, 
< >j7j= 0. Therefore, for any given point x G M, one can choose a normal coordinate system 
such that the following relation holds: 

y^djf{x,y) = (4.4) 

at that point. The existence of normal coordinates is of fundamental importance in the sequel. 
For instance, one can get an expression for the covariant derivative of V along the integral curve 
of V, 

< ^D>yV = {V^d,V')-^, (4.5) 

using a coordinate frame {g^, k = 0, ...,n — 1}. Note that this expression is not a partial 
differential equation because it holds only at the point x. 

Developing the above expression further, we have the relation: 

< >y Vix) = ^ I dvolix,y)y^f{x,yy 

V0l{supp{fx)) JsuppiU) 



■dj{ ^ / , dvol{x,y)y^f{x,y)). 

voUsuvpifT.)) Jsuvvif^) 



' V0l{supp{fx)) . Jsupp{fx) 

supp{fx) is compact and small compared with E{x), for each x G M. It is the right hand of this 
equation which we should estimate. 



dvol{x,y)y^fix,y)dj{ ^ [ _ dvol{x,y)y^ f {x,y))\ = 



vol{supp{fx)) ^JsuppiU) vol{supp{fx)) . Jsupp{fx) 



dvol{x,y)y^ f{x,y) ( - 



vol{supp{fx)) ^ JsuppCu) voP{supp{fx)) 
■dj{[ dvol{x,y)f{x,y))))+——^—Y—([ dvol{x,y) y^ f{x,y) ■ 

Jsuppih) ' V0l{supp{fx)) ^ JsuppCh) 

■dj{ [ , dvol{x,y)ff{x,y))) = 

Jsupp{fx) 

— — y— ( / dvol{x,y)y^f{x,y)dj{ [ dvol{x,y) f{x,y) ) </>) + 

voV{supp{fx)) ^ JsuppiU) Jsuppif^) ' 

dvol{x,y)y^ f{x,y)dj { _ dvol(x,y) y^ f{x,y)) ). 



VOp{sUpp{fx)) ^ Jsupp{fx) Jsupp{fx) 
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Shifting the variable of integration — < y'^ >= —6^{x,y), one obtains for the above expres- 
sion 



< >y V{X) 



1 



VOp{sUpp{fx)) ^ 



dvol{x,y)y^ f{x,y)- 



d 



i su 



■jy f dvolyf{x,y) )</>) + 

SUpp(fx) 



^ ( / dvol{x, y) y^ f{x, y) dj ( / dvol{x, y) (< y^ > +S''{x, y)) f{x, y) 



Jsupp{fx) 

,9/~~7TTt(/ dvol{x,y)y^ f{x,y)dj{ I dvol{x,y)5^{x,y) f{x,y))\ 

voV{supp{fx)) ^ JsuppUx) Jsum(U) ' 



Since y^dif{x,y) = at the origin x of a normal coordinate system and since fx is a smooth 
function of y, we can Taylor expanded the above expression: 

— — ^— y— ( / dvol{x,y)y^ f{x,y)dj{ [ dvol{x,y) S''{x,y) {f{x,y)+^ {y^-y^ 

Since the averaged connection is an affine connection, there is a coordinate system such that 
< >jfc= 0. This is reflected in the expression for the averaged Vlasov equation, which has 

the form y^ djf{x,y) = at one given point x G M. Then, we get the following for the above 
expression 



( < >^ V{x)) 



VOp{sUpp{fx)) ^ JsuppiU) 

d,{[ dvol{x,y)5\x,y)%{y' -y'))). 



dvol{x,y)y^ f{x,y) 



'supp{fx) 

{y^ ~ y') and 6''{x, y) are bounded by the diameter a (remember that in taking the moments we 
have substituted f^, supp{fx) by {fx, supp{fx)). Therefore, 



SUppifx) 



dvol{x,y)y^ f{x,y) ■ 



< 



VOp{sUpp{fx)) 

■dj ( / dvol{x, y) 5\x, y) % {y' - y')dk ) ) 

JsuppUx) ^ 

( / dvol{x,y)y^ f{x,y)dj ^) 

^ Jsupp(fx) ^ 

dvol{x, y) 5^{x, y) (y' - y^)dk ) + 

/ r) 

( / _ dvol{x,y)y^ f{x,y)^) 
^ Jsupp(.L) 9yt ' 



voP{supp{fx)) 



SUpp{fx) 
1 



VOp{sUpp{fx)) 



< 

1 



SUpp{fx) 
1 



dvol{x,y) dj5^{x,y) (y' - y')^^) 



voP{supp{fx)) 



< 



df 



SUpp{fx) 



dvol{x,y)y^ f{x,y)dj 
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(/ dvol{x,y)\\5\x,y){y'-y')dk ) + 

^ dvol{x,y)y^ f{x,y) 



SUpp(Jx) 



VOp{sUpp{fx)) 

■if _ dvol{x,y)dj\\5\x,y){y'-y')dk\\-) 



'supp{fx) 

One can find a bound for each of these integrals. For instance, using the Hoelder inequality for 
integrals [13] 



U<t>dy.\<[j Wf\f Wf\ - + - = 1, l<p,q<oo. 



!/<? 1 1 
- + - 

/X ^ ^ P Q 

We will use this inequality several times for the case p = q = 2. Then, we obtain the following: 



SUppifx) 



< 



dvol{x,y) {y'-y')5\x,y)dk )<( [ dvol{x,y) \{y'-y')\ \\ 5\x,y)dk \\-) 

dvolix,y)\{y'-y')fY- ( / dvol{x,y)\\S\x,y)dk\\iy < 

< vol^{supp{fx)) a ■ ( / dvol{x,y) 11 S''{x,y)dk 11^) ^ 



< 



SUpp{fx} 



The Euclidean volume was 

volE{supp{fx)) ■■= / _ dvol{x,y). 

J SUpp{fx) 

To bound the second factor we use the following [6] : 

1 3 

\\S{x,y)\\f^ <\\<y>{x) -y\\f^ < \\e + y-y\\f^ < \\e\\fi + \\y - yWa < +a = -a. 

The parameter e is related to the norm of the mean velocity and y is in the support of the 
distribution /. Therefore, a bound on the integral is 



dvol{x,y) {y^ -y^)5^{x,y)dk^ 



'SUpp{fx) 

Similarly one obtains the following bound: 



< 77 • volE{supp{fx)) ■ a 

ri I 



([ dvol{x,y)dj{8\x,y)){y'-y')dk) _< (f _ dvol{x,y)\{y' - y'fY ■ 

dvol{x,y)^dj{5^{x,y))dk\^.Y < 
< vol^{supp{f:a)) a ■ ( _ dvol{x,y) \\ dj{6'' {x,y))dk\\^-Y . 



supp{fx) 



Isuppifx) 

Because of the definition of the corresponding Sobolev norm one obtains 



dvolix,y) \\djiS'{x,y)) {y'-y')dk\\-) < voll{supp{U)) a ■ ||5,<5^||o,2 



SUpp(fx) 
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= voll{supp{h))a-Y,¥%log{5l)\\o,2< voll{suppU^^^ 

k k 

Similarly, 

( / _ dvol{x,y)y^ f{x,y)dj —:t)\ < y2 { / - /(x,y)|^ ) 

.(/ «(.,,)|a,£?fjl^>r)*. 

^ J supp(fx) ' 

The second factor is equal to the Sobolev norm ||i9j/a;||i,i- The first factor is bounded in the 
following way: 

( / dvol(x,y)\y> j{x,y)'^y <( dvol{x,y)\y>\^ f{x,y)y = 

= vol{supp{f)) {< >)i 

Therefore, we get the bound: 

_ dvol{x,y)y^f{x,y)dj-^)\ < vol{supp{fl)) [Y.(< ly^'' >)'^ ■ H^,- ) 

In a local frame where the vector field < U >= (?7°, 0), this contraction can be re- written as 

(< {yy >)^ • l|5o/.||i,i = ||(< |yT >)^ • doMi,i = ll(< iy'? >)^ • dofAi,i = 

= ||(< {y^f . {dohf >)^||i,i = ||(< {yi . djflf >)^||i,i. 

Using normal coordinates system associated with the affinc connection < > and using the 
averaged Vlasov equation, one obtains that ||(< {y^ ■ djfxf' >)2 ||i,i = 0. 

Finally, we can bound the following integral 



SUpp{fx 

n-1 



dvol{x,y)y^ f{x,y) —j) <( dvol{x,y)\y^ f{x,y) < 

dy ^ ^ Jsuppifx) ^2/ ^ 



<Y.i[ - dvol{x,y)Wf{x,y))\^y-{[ _ dvol{x,y) l^^^l'^Y . 

„_n J sunp(fx) ^ Jsuppifx) 



j=0 supp{fx) ^ Jsuppifx) 

Like in the previous integral, we get: 



if [ dvol{x,y)y^ f{x,y) ^)\ < vol{supp{fx)) ■ (< \y^\'^ >)2 • ||/^| 

^ Jsupp( fx) oy' J 



1,1- 



ISUpp{fx) 

Using these bounds, we obtain the following relation: 



<LD>yV{x)\U<——^—j- (! _ dvol{x,y)y^ ~f{x,y)d^^) 

voV{supp[fx)) ^ Jsuppifx) "y ^ 



supp{fx) 



3 

■-V0lE{supp{fx))0? + 
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n—1 „ J, 

df 



TTTT ( / - dvol{x,y)y^ f{x,y)—j 

■voll{supp{U))-a^-{Y. ll^iM'5')llo,2) < 



k 

n—1 



- 7YV^ ■ vol{supp{f^)) ■ ( V(< W? >)^ • ||/x||i,i • vol%{supp{f:c))- 

VoV{supp[fa:)) 



k 

n-1 



j=0 k 

In a local frame where the vector field U has components (?7°,0), one has 

n-l 

E« \y'\' (E \\93lo9{5l) ||o,2 = (< |yT >)^ (E 11^0^05(5^) ||o,2 = 

3=0 k k 

k 

Note that the normal coordinate system (that we are using) coincides with the adapted coor- 
dinate system, associated to the vector field U = (?7°,0. In this coordinate system, there is a 
bound < y'' >< 1 + a. Therefore, 

n-l 

E « l^'l' »^ (E ll^iM^') Il0,2 < (< >)^ • (E 11^0 M-^^) Il0,2 • (1 + «). 

3=0 k k 

We get the following result, 



II < >y Vix)\\, < ^^^Eisuppih)) 115^^,^(5.) 11^^^) . 11/^11,^, ^o{a% 

vol{supp{Jx)) ^ 

□ 

Usually we are interested in the behavior of the fluid in compact regions of the space-time M: 



Corollary 4.4 For compact domains K C M and under the same hypotheses as in theorem 
4-3, the following relation holds: 

II < >y V{x)\\f^ < -CiK) ■ + O(a^), 

for some constant C(K). 

Proof: Take for the constant C(K) to be 

1 

C7(K) = max^^A (E Wo9{5l)\\,,,) ■ \\h\\,,, ]. 

i- vol{supp{fx)) ^ J 

□ 

These bounds axe quite general are meaningful in the case that 1 >> a. 

In the preceding results the ultra-relativistic limit {E » 1) was not essential, although it 
in intrinsically used when we write conditions on the smoothness conditions of supp{f) and the 
measure used. 
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Let us consider the normalized averaged velocity vector field: 

V 



u 



Since < > docs not preserve the Minkowski metric 77, the covariant derivative of u in the 
direction of u using the Lorentz connection is 

<'^D>^u= -^J^ <Ljj>^v+^(v- {log{ri{V, V))))v. (4.6) 

The first terms is bounded by theorem 4-3, since r){V,V) > 1. The total derivative of r){V,V) 
along a trajectory of V is 

V)) = jC.y{ri{V, V))=V- (r?(y, V)) = 2r?( < ^L>>^ f) + ( < >y r?) {V, V). 

As we have proved the first term is of order a^. Using normal coordinates for < >, one can 
compute the second term, 

( < >yri){V,V) = ri{V,V)F^rrr<S"'{x,y)5\x,y)d^{x,y) > V^'VsVi. 

We can estimate these contributions: 

Proposition 4.5 Under the same assumptions as in theorem 4-3, the following relation holds: 

1 

<^D>,u< ^"f'^PPf;^^ \\dolog{S'.)\\o,2) -Wflhi-a^ + 0{a% (4.7) 
vol{supp{fx)) ^ 

where n is the dimension ofM. 

Proof: The first term of the right hand of the equation 4-6 is bounded by theorem 4-3. The sec- 
ond term is bounded using Hoelder's inequality for integrals in the following (this is a particular 
case of [16, pg 62]): 



/ dvol{z) h{z) ■ ■ ■ fm{z)\ <n ( / dvol{z)\h{z)\'"'y\ ^Pfc = l, I < Pk < 

^ S 1. 1 *J S 7- 



00. 



fc=l 



In particular one can apply this inequality to the third order moment < 5"^ (x, y) 5" (x, y) S {x, y) >: 

<S^ix,y)5%x,y)5\x,y)> = ^ [ dvol{x,y) f{x,y) S^{x,y) S%x,y) 5\x,y) 

volisuppifx)) Jsuvv(L) 



< „ ^ ,,,, ( / dvol{x,y)\f{x,y)S"'{x,y)fy-( [ dvol{x,y) \f{x,y)d%x,y)\'y ■ 

V0l{supp{fx)) ^Jsuppif^) ' ^Jsuppif^) ' 

■{ inf dvol{x,y)\j{x,y)b\x,y)'fY . 

Since the distribution function is positive on supp{f) and fx < 1. Then, since equation (3.1), 
one can substitute in the integrations fx — > fx, which implies that 

<6^{x,y)6%x,y)6'{x,y)>\ =0{a^). 

Since the norm fj{< y >,< y >) > 1, one gets a third degree monomial term in a for the 
covariant derivative < >{. u. □ 
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Corollary 4.6 Under the same assumptions as in theorem 4-3, in a compact domain K C M 
one obtains the following asymptotic relation: 

< ^D>uu< C(K) ■ o? + 0{o?), (4.8) 

/or a convenient constant C(K). 

Let us consider a local Lorentz congruence, which is a set of auto-parallel curves of the Lorentz 
connection ^D, for a set of initial conditions at each (tQ,x), x € M^q where Mf(, ^ M is a 3- 
dimensional spatial sub-manifold. One can consider in a similar way the congruence associated 
with the averaged Lorentz connection with the same initial conditions. Note that, while the 
Lorentz connection preserves the Lorentz norm r]{x,x) of the tangent vectors of the geodesies, 
this is not the case for the averaged Lorentz connection). 

Theorem 4.7 Let F be a closed 2-form and the associated non-linear Lorentz connection. 
Then the solutions of the Lorentz force equation ^VxX = (i^F)'' can he approximated by the 
integral curves of the normalized mean velocity field u of the distribution function f{x,y), where 
f{x, y) is a solution of the associated Vlasov equation ^xf = 0- difference is controlled by 
polynomial functions at least of order 2 in a: 

\\^Duu\\ < a2{x)a^ +0{a^). (4.9) 

Proof: We repeat an argument that we have used before. By proposition 3.1, both distribution 
functions / and /, solutions of ^xf = ^^'^ < > ^i-^e such that: 

\f{t,x{t),xit))-fit,xit),xm < {Cix)\\F\\ix)Cl{x){l + B2ix)a))a''E-H^ + 

+ {k{x)\\F\\{x) kI{1 + D2{x)a))a'^ E-^t. 

Therefore, the corresponding mean velocity fields are nearby as well, because the linearity of the 
averaging operation and because of the relation: 

\f{t,x{t),x{t))\ - \f{t,x{t),x{t))\ <\f{t,x{t),x{t))-f{t,x{t),x{t))\ < 

< {C{x) ||F|| (a;)C| (x) (1 + B2{x)a)) E'^ t^ + {k{x)\\F\\{x) K| (1 + D2{x)a)) o? E'^ t. 

Then, their corresponding integral curves and the associated local congruences are also numeri- 
cally similar. By corollary 4-6, for narrow distributions, the normalized mean field u associated 
with /, 

II < >s u\\ < 0.20"^ + O(a^). 

for some known function 02 (x). Remember that, one can interpolate smoothly between the 
connections and < >. Therefore, locally, one can interpolate smoothly between their 
integral curves. Also, because the smoothness of the solutions of the geodesic equations in the 
parameter of interpolation, there is a function a in a small open neighborhood of M such that 

^DuU < a2{x)a^ +0{a^). 

□ 
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5 Discussion 



Theorem 4-5 and Corollary 4-6 show when the charged cold fluid model in the description of 
the dynamics of a collection of particles interacting with an external electromagnetic fleld in the 
ultra-relativistic regime. It is particularly interesting that we have proved the result without 
using additional hypothesis on the behavior of the higher moments of the distribution, except 
that the distribution is narrow and some convenient smoothness conditions. 

There are some technical issues that we would like to mention. 

1. The distribution functions / and / are at least C^. However, let us consider a Dirac delta 
distribution with support invariant by the flow of the Lorentz force [6]: 

f{x,y) = ^ix)Siy-V{x)). (5.1) 

Since the width of the distribution is zero, a = 0. One can use this distribution as a 
solution of the Vlasov equation in two space-time dimension spaces, for a proper value 

of the function ^ . This example and the fact that the bounds found in section 4 arc 
formulated using Sobolev norms suggests the possibility to generalize the results to bigger 
spaces. 

2. The need for the condition supp{fx) C supp{fx) needs to be investigated. We have seen 
that one can reduce supp(fx) conveniently. Indeed, the inclusion is non-trivial, since the 
delta function limit is a solution of both, the averaged Vlasov equation and the original 
Vlasov equation. On the other hand, the integro differential equation < > / = is 
in some sense simpler than the original Vlasov equation ^x(/) = 0- This is because the 
averaged connection is simpler than the Lorentz connection. This fact provides an insight 
why supp{fx) C supp{fx). 

3. The same method can be applied to other fluid equations. Therefore, depending on the 
specific bounds, one can decide which model is better in a particular situation 

4. Although the results presented in this work are Lorentz covariant (in the sense that the 
normal coordinate system has an invariant definition) , they are written in a non-covariant 
way. This is because wc have introduced norms associated with particular reference frames. 
These norms arc essential in our approach. At the moment, wc don't know how to give a 
Lorentz invariant presentation of the content. Therefore, further research is necessary on 
this direction. 
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